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ABSTRACT 


The  unsteady  lifting  surface  approach  is  utilized 
for  the  marine  propeller  case  and  the  corresponding  surface 
integral  equation  is  solved  for  the  Weissinger  mathematical 
model.  The  applicability  of  the  Weissinger  method  to  the 
nonstationary  flow  case  is  studied;  The  kernel  function  is 
expressed  in  closed  forms  after  some  mathematical  simplifica¬ 
tion.  Prom  numerical  calculations  of  propeller  loading 
which  are  restricted  to  a.  four-bladed  propeller  of  sector 
type  blade  form  with  different  blade- area  ratios  and  various 
pitch- diameter  ratios,,  conclusions  are  drawn  as  to  three- 
dimensional  effects  as  well  as  to  the  dependence  of  propeller 
loading  on  the  important  parameters,  as  the  blade-area  and 
pitch-diameter  ratios. 
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INTRODUCTION 


The  ship  hull  induces  a  disturbance  flow  upon  the 
propeller  plane.  Although,  in  general,  this  flow  is  time- 
independent,  its  distribution  on  the  propeller  plane  is  non- 
uniform,  and  hence  the  rotating  blades,  which  cut  through 
this  disturbance,  experience  time-dependent  gust.  If  this 
disturbance  flow  is  expanded  into  angular  harmonics,  each 
component  exerts  a  sinusoidal  gust  upon  the  blades.  Using 
Sears'  two-dimensional  response  function  for  the  sinusoidal 
gust1  in  a  stripwise  manner,  Breslin  et  al2  evaluated  the 
oscillatory  thrust  and  torque  of  the  propeller  due  to  this 
kind  of  disturbance.  However,  correlation  with  existing 
experimental  results  was  very  poor,3  and  they  attributed  this 
to  application  of  the  two-dimensional  method  to  a  typically 
three-dimensional  problem. 

At  this  stage,  it  seems  important  to  try  a  three- 
dimensional  approach.  The  three-dimensional  lifting  surface 
integral  equation  for  the  screw  propeller  in  a  steady  flow 
was  derived  by  Sparenberg.4  A  slight  modification  gives  the 
equation  for  the  non-steady  state.  The  nature  of  this  equa¬ 
tion  was  examined  by  Hanaoka.5 

Owing  to  the  complicated  form  of  the  kernel,  how¬ 
ever,  It  seems  difficult  to  treat  this  equation  in  an  exact 
form.  In  this  report,  by  adopting  the  Weissinger  approxima¬ 
tion,  the  surface  integral  -equation  is  converted  into  a  span- 
wise  line  equation.  The  applicability  of  the  Weissinger 
method  to  the  three-dimensional  propeller  problem  in  the 
steady  condition  has  been  well  examined,  but  applicability 
to  the  non-steady  problem  remains  to  be  tested.  In  Appendix  A 
the  two-dimensional  sinusoidal  gust  problem  is  treated  under 
the  Weissinger  approximation  and  the  result  is  compared  with 
Sears'  exact  value.  The  result  is  also  applied  stripwise  to 
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the  present  problem  and  a.  comparison  with  the  results  of  the 
three-dimensional  Weissinger  approach  is  shown.  This  com¬ 
parison  (both  approaches  are  obtained  under  the  same  Weissinger 
approximation)  will  provide  important  knowledge  on  the  three- 
dimensional  effect  in  the  non-steady  state. 


FORMULATION  OF  PROBLEM 


In  this  section  the  derivation  of  the  integral 
equation  for  the  .angular  harmonic  component  of  the  disturbance 
flow  is  shown.  Use  is  made  of  the  acceleration  potential 
method.  The  necessary  notations  concerning  the  coordinate 
system  are  shown  in  Fig.  1. 

The  pressure  field  Pp(x,r,.0)  due  to  a  unit  pressure 
pole  1 is  given  by 


P  +  r2  +  p2  -  2r  p  cos  ( 9-<f> ) 

If  the  pole  is  pulsating  in  its  strength  with  frequency  cd  and 
rotating  in  the  negative  direction  of  9  with  angular  velocity  £ls. 

i  icot 

p  _  1  — ■ - § — - — . - . - - - - - -  ^  2^ 

p  ^  ^(x-l)2  +  r2  +  p2  -  2r  p cos 

The  linearized  relation  between  the  perturbation  velocity 
potential  ®(x,r,0,t)  and  the  corresponding  pressure  field 
P(x,r,^,t)  in  the  main  (x-.direction)  uniform  flow  U  Is  given  by 


U3x  Bt 


(3) 


where  /j  is  the  density  of  the  fluid.  If  it  is  assumed  that 
0  =  0  at  x  =  -«>,  the  solution  <t>  for  given  P  is 


x 

P(r’,r,0,t  + 
-  00 


dr’ 


(4) 
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On  substituting  eq.  2  into  eq.  4  and  putting  r’  =  t+|,  the 
velocity  potential  corresponding  to  of  eq,  2  becomes 


lost  r 


x- 


4>  = 

P 


e _ 

47TpU 


ig(T-X+t) 


+ra+p2-2rpcos^@o-nt-^(T-x+^) 


■dr 


(5) 


When  the  observation  point  (x,r,$)  is  rotating  together  with 
the  pressure  pole,  putting  <£  =  0Q  -  fit. 


4>  =  - 


lost 

e 

4TTp.tr 


i^j(r-x+£) 


'\Jtz+vz+pz  -2rpcos  j~  (  eo-0o)  -fi  (  t-x+£ 


U 


>] 


dr 


(6) 


If  the  pressure  pole  is  replaced  by  a  dipole,  the  correspond¬ 
ing  velocity  potential  <&d  is 


*d  " 


ax 


(7) 


where  - — -is  a  differential  operator  with  respect  to  (£,  p,  0  ) 
^X 

In  the  direction  of  the  dipole  axis.  Velocity  induced  by  this 
dipole  is 


Wd  (x,r,0Q) 


_a_ 

^X 


4>  = 

»  d 


ax’  ax 


p 


(8) 


where- — •  is  an  operator  with  respect  to  (x,r,$  )  in  the 

ax'  0 

designated  direction  of  the  velocity. 


In  the  acceleration  potential  method,  the  lifting 
surface  is  regarded  as  a  mono-layer  of  pressure  dipoles,  the 
axis  of  which  is  normal  to  the  surface  and  the  strength  of 
which  is  proportional  to  the  pressure  jump  across  the  surface. 
Therefore,  it  is  easily  seen  that  eq.  8  is  the  elementary  com¬ 
ponent  of  the  lifting  surface  integral  equation  for  the  simple 
harmonically  time-dependent  problem  of  a  rotating  system. 
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The  phenomena  appearing  at  corresponding  points  on 
all  blades  due  to  the  qth  harmonics  of  the  disturbance  flow 
are'  expressed  in  the  following  simple  form: 

tn  -  C0  e1^  -  (9) 

where  =  -^(  n  -  l) ,  n  is,  the  name  of  the  blade  n  =  1,  2, 

. . . . .N,  and  N  is  the  total  number  of  blades .  Summing  up  the 
effect  of  all  blades  upon  a  particular  blade  (n  =  l) ,  and 
considering  that  the  geometrical  equation  of  the  lifting  sur¬ 
face  of  the  screw  propeller  is  expressed  by  • 


=  |  ( or' i  =  i  0 J  where  a  =  #, 


a  o' 


Q 

U' 


(10) 


leads  to  the  integral  equation  for.  the  present  problem  which 
can  be  written  as 


iq(flt  -  <(>  ) 

V(  r)  e 


iqftt 


=  e 


K(r,  (f>  ,  p,  9  )  d0  dp  (ll) 


Here  S(  p,  0  )  is  the  lift  distribution  on  the  blade,  and 
K(r,$o,p,@o)  is  the  kernel,  which  is  obtained  in  the  following 
form  by  using  eqs .  6  to  10: 


N 

£  lim 


K(riV,ao)  -  - 

6  =  !('Veo)-(x-^-+0 


_c> _ a_ 


iq[-en+a(T-x+t)] 


—  —  —  —  —  —  —  —  =-dT 

/T2+r2+p2-2rpcos  [(0o-<4o)+en  -a(T-x+£)] 


(12) 


The  directions  of  the  differentiations  in  —  and  —  are 

&X  ^x  * 

perpendicular  to  the  lifting  surface,  that  is,  from  eq.  10 
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a 

r 

L  S 

i 

a 

dX* 

^ l+a2r2 

\ 

r2 

a  $ 

d 

P 

/„  a 

1 

a 

Sx 

\]l+a2  p2 

\  H 

P2 

ae 

The  limiting  process 

lim  6  =  (x-£)  )— >0 

\  '  T-  o  O 


(13) 


is  introduced  in  order  .to  avoid  the  difficulty  due  to  the 
singularity  in  the  mathematical  manipulations.  Physically* 
it  means  that  the  manipulations  are  performed  taking  the  con¬ 
trolling  point  on  a  surface  slightly  shifted  from  the  lifting 
surface  and  that  finally  the  former  surface  is  brought  in¬ 
finitesimally  close  to  the  latter. 


SOLUTION  BY  THE  WEISSINGER  APPROXIMATION 

APPLICATION  OP  THE  WEISSINC-ER  METHOD 
AND  SOLVING  PROCEDURE 

In  the  Weissinger  method*  the  surface  distribution 
of  the  lift  is  replaced  by  a  line  distribution  along  the 
l/4-chord  line*  and  the  controlling  points  are  taken  along 
the  3/4- chord. line .  This  method  gives  an  exact  lift  for  the 
two-dimensional  stationary  wing  problem,  and  it  is  applicable 
also  to  the  stationary  three-dimensional  wing  problem  under 
the  assumption  .of  not  too  low  aspect  ratio.  In  fact,  in  the 
stationary  propeller  problem,  it  has  been  shown  that  this 
method  can  describe  the  effect  of  the  helical  free  vortex 
fairly  well . 

However,  in  the  nonstationary  case,  there  is  another 
kind  of  free  vortex  which  is  parallel  to  the  span  ( in  the 
propeller  problem  the  radial  free  vortex) ,  and  it  seems  not 
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to  have  been  clarified  how  precisely  the  Weissinger  model  can 
obtain' the  effect  of  the  latter  kind  of  free  vortex  The 
applicability  of  this  method  to  the  nonstationary  problem  is 
tested  in  Appendix  A  by  treating  the  two-dimensional  sinus¬ 
oidal  gust  problem.  The  result  indicates  that  this  method 
can  be  used  up  to  a  reduced  frequency  of  1.3*  which  corre¬ 
sponds  to  blade  frequency  harmonics  (a  =  N)  in  a  propeller  of 
area  ratio  0.4. 


This  approximation  is  based  upon  the  first  and 
second  terms  of  Birnbaum’s  chordwise  lift  distribution*  and 
hence  on  the  left-hand  side  of  the  integral  equation  the  cor¬ 
responding  components  should  be  picked  up.  In  the  two- 
dimensional  stationary  case*  this  can  be  done  exactly  by 
applying  Glauert's  lift  operator* 

3  ( 1  -  cos  a) da  where  x  =  -b  cos  a*  x  is  the  chordwise 
o 

coordinate  and  b  is  the  semi-chord  length.  The  same  operator 
is  tried  for  the  two-dimensional  but  nonstationary  case  in 
Appendix  A*  and  the  result  indicates  its  satisfactory  applica¬ 
tion  to  the  nonstationary  case  up  to  the  reduced  frequency 
1.3*  In  the  following,  treatment*  the  same  procedure  will  be 
used  for  the  nonstationary  three-dimensional  case  under  the 
assumption  of  a  not  too  low  aspect  ratio. 


Thus*  the  application  of  the  Weissinger  method  con¬ 
verts  eq„  11  (the  surface  integral  equation)  into  the  follow¬ 
ing  line  integral  equation: 


W(  r) 

where 

W(  r) 
L(p) 


=  )  L(p)K(r*p)dp  * 

0 

-7T  - 

l  iqe^cos  a 

-  V  (r)  J  e  (l  -  cos  a)  da  * 

0 

is  spanwise  line  distribution  of  lift* 


(14) 

(15) 


K(r*p)  =  [K(r*<jf>o*p*0o)J( 


-  9  =  0, 


(16) 
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and 


0^  is  semi-chord  length. 

The  integral  equation  (eq.  1*0  can  he.  treated  in 
the  following  way.  By  dividing  the  span,  into  intervals  of 
equal  length  21,  and  regarding  the  lift  distribution  on  each 
fraction  as  uniform.,  eq.  16  can  be  written  in  the  form  of  a 
set  of  simultaneous  linear  algebraic  equations 


(17) 


The  element  of  the  kernel  matrix  is  given  by 


rV^ 

Kpq  =\  K(rP  •  p)dp  •  (l8) 

K  -i 

q 

The  singularity  of  eq.  12  is  involved  in  the  element  K  ,  but 
that  singularity  is  integrable  in  the  p  integration.  Thus,  the 
kernel  matrix  element's  are  free  from  singular  behavior,  and 
the  spanwise  lift  distribution  for  given  W^,  and,  in 
addition,  the  inverse  matrix  °^[Kpq]->  can  obtained  without 
any  difficulty.  In  the  nonstationary  problem,  all  the  ele¬ 
ments  in  eq.  17  are  complex  numbers.  Hence,  if  the  span  is 
divided  into  P  parts,  the  kernel  matrix  is  of  the  order 
2P  x  2P,  but  a  conventional  digital  computer  can  easily  give 
the  required  results. 


SIMPLIFICATION  OF  THE  KERNEL 

By  putting  y  =  -a (r  -  x  +  |) ,  interchanging  the 
order  of  the  integral  sign  and  differential  operation  in 
eq.  12  becomes  possible.  The  term  within  the  brackets  of 
eq.  12,  designated  now  by  H,  becomes  on  substituting  eq.  13, 
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00 


H 


1 

a 


-iq( en+y) 


d  1_  J3 

/lAA"2  '  ~dx  r2  b<p'o  yi+aV'  ~&€  ■  /S0O 


V[|y"^x~^)]2  +  1,2  +  p2  -  Srpcos  [(eo-,0o)  +  0n  +  y  J 


rdy  (19) 


Further,  If 


H 


-iq( 0n+y) 


-y  -  (x-g)  =  X  and  @n  +  ( 0Q  -  0Q)  +  y  =  0  ,  ( 20) 

~-r  (a  A  -  A  A) _ E _ (a  A  -A  A 

V l+a2r2  5x  r2  50  V 1+a2  p2  P2 

dy 


V  X2+r2+pa-2rpcos© 


=  +  (e  -  +  y 


n 


X  =  |y  -  (x-|) 


Changing  the  Integral  variable  from  y  to  0 


H  =  "l 


-iq[0  -  (0O-0O)] 


0  +(0  -0  ) 
n  o  o 


_ /ad  1  5 


/ad  1  5  v 

'  *N  O  V  *  ' 


V  1+a2  r2  dx  r2  d0  yi+a2p2  dx  p2  d0 


where 


y  X2+r2+p2-2rpcos© 
1 


x  -  i  9  -  |  5n  +  6 


d© 


a.  "  (x  ~  €) 


(21) 


(22) 


Here,  under  the  assumption  that  the  pitch  of  the  propeller  is 
low — i.e.,  a  Is  large — the  following  two  approximations  are 
Introduced: 
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(23) 


1)  (a  JL  „  JL_ 

nj l+aar2  dx  r2  .be  bx 
- £ — ~  (  „  S  1  S  \  .  5 

V1+a  ap2  pa  Se  Sx 

2)  the  change  shown  in  Fig.  2  of  the  integral  path  in 
eq.  22  from  the  original  one'  expressed  by  light  solid 
lines  to  the  heavy  solid  lines. 

These  two  approximations  give  K,  eq'.  12,  in  the  following  form: 
1  "lq®b 

K  =  Wju.Ua  e  [Kx  +  Kj  , 
where. 


N 

j 

■§n+(  V2> 

b 2 

e-iq© 

n=l 

lim 
6->0  . 

-(0  -e  ) 
n  ^0  o' 

L 

Vxa+r2+p2-2rp 

cos  0 

K2=  lim 
6->0 

00 

£ 

m=l 

r 

e-lqe  1 

Jo 

bxz 

)Jxz+v2+p2~2vp  cos  0j 

X=mi-0 
a  c 

(24) 

S  _  27T 
o  N 

and  ($O_0O)  =  in  the  Weissinger  approximation. 

The  approximations  l)  and  2)  are  rather  rough  ones 
if  they  are  adopted  separately,  but  their  combined  use  is 
justified  by  a  clear  geometrical  meaning  as  well  as  good 
numerical  results.  This  point  is  discussed  in  Appendix  B. 

It  is  also  worthy  of  notice  that  K1  depends  upon  0^, 
the  semi- chord  length,  but  is  independent  of  a,  the  inverse 
of  pitch.  On  the  other  hand,  K2  depends  upon  the  pitch  but 
is  independent  of  the  chord  length.  This  fact  is  a  great  con¬ 
venience  for  the  numerical  calculations. 
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EVALUATION  OF  THE  KERNEL 

Kx  .  When  p  £  r,  the  limiting  process  <5->0  can  be  taken 
before  the  integration,,  then 


Kj.  =  £ 
n=l 


N  rj„+<V2> 


-iq® 


e  -0. 

n  b 


(  ^r2+p2-2rp  cos  ®)£ 


d® 


(25) 


Further, 


s-iq0 


Q  _Q  yr2+p2-2rp  cos 
n  b 


•d® 


=  _JL_ 

(r+p)3 


-i/2(Tr+en-ey2) 


1  -i/2(  7T+en-0b) 


,i2q^ 


^/l-k2sin2^ 


(26) 


where. 


f  =  -  i  ( 71+®)  and  k2  = 

2  (r+p)  = 


(27) 


As  Is  shown  in  Appendix  C,  the  integral 


u 


i2q^ 


Vs 


^1  -  k2sin2^ 


■d^ 


(28) 


can  be  evaluated  by  the  following  recurrence  formulas 


C  = 

Ln+1  2n-l 


(1  -  4-)  I  -  (USt)  I 


k2 


n  '2n-l  n-1 


-l 


2  2 


=i2n^ 


2n-l  ,2 


k2  y1  l-k2sin2^ 
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and 


Zo  = 


1  -  kc 


E(  Kf) 


k2  sin2^ 

T 


yi-k2sin2^/  J 


u 


f 


l 


I  +—  P(k^) 
°  k2 


t 


u 


+i 


_2_  1 
k2  ^  l-k2sin2^/ 


t 


(29) 


where  F(k,$)  and  E(k,^)  are  elliptic  integrals  of  the  first 
and  second  kind,  respectively.  In  the  case  of  p  /  q,  the 
integration  with  respect  to  p  in  eq.  18  can  he  evaluated  by 
conventional  numerical  methods  or  by  taking  some  conventional 
mean  value . 


In  the  narrow  band  expressed  by  p  -  r  <<r  ,  which 

ir  I  Jr 

exists  in  the  element  p  =  q,  a  special  treatment  is  necessary. 
As  is  shown  in  Fig.  3*  this  area  is  divided  into  the  following 
three  regions: 

1)  r  -  e  <  p  <  r  +  e  ((Xe<j31  02)  and-y<0<y  (  0<ry«l) , 

p  o  p  o  i 

which  include  the  singular  point  p  =  r  0=0,  (double 
hatched  region)  p 

2)  Tp  -  eQ<  p<  r^  +  eQ  but  excluding  -  y<Q<y,  (dotted 
region) 


3)  rp  -  /3i  <  p  <rp  -  eo  and  rp  +  eQ<  p  <  rp  +  /32,  (  single 
hatched  region) 


Designating  the  contribution  of  these  three  regions  of  K1  as 
Jlf  J2  and  J3  respectively,  eq.  24  gives 


Jj_  =  lim 

6=X->0 


•rp  +  £c 


>r  -  e 

P  o 


,  -Y 


-7 


•iq© 


cb?2  \/x2+r2+p2-2rpcos0 


■  d0  dp  (30) 


R-940 

-11- 


Jo  = 


lim 

5=X-»0 


r  +  e 
P  o 


e  N  re  +1/20 

o  ^  £  n  o 


+ 


“6b  4-7 


n_1  Je  -0, 

n  b 


dx2  Vx2+ra+p2-2rpcos© 


-d0  dp 


(31) 


J3  =  lim 
X=6->0 


F  r  -  e 

p  o 


4- 


L  “rP  -  ^ 


r  +  fiP 
P  Pa 


r  +  £ 

P  o 


r  A +  Vae0 
£ 


C 


-iq© 


0-0, 
n  b 


dx2  y  X2+r2+pa-2rpcos0 


d©  dp 


(32) 


The  most  convenient  form  for  evaluating  the  total  contribution 
J  of  these  three  parts  is  the  limit  value  expression: 


J  =  lim  (J  +  J  +  J  ) 

e— >0 
o 


(33) 


As  is  seen  in  Appendix  D: 


lim  J  =  lim 

e— *-0  e — >0 

o  o 


r  e 
P  o 


(34) 


The  integrand  of  J2  is  finite  through  the  whole  integration 
range ^  and  hence 


llmn  Ja  =  0 
£—>0  2 
O 


(35) 


In  Appendix  E,  J  is  given  in  the  following  form: 
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where  0  and  /32  should  be  taken  to  satisfy  the  relation 


k{  = 


where 


1 1  t 

q 


2rP“/3i  2ry^2 
-1/2(7 r-ej 


-T7 


-1/2(77+0^2) 


(37) 


N  f -1/2(77+0  +0  /2)^ 

+  £  \ 

J_1/2(7T+0n-0b) 


n= 


y  l-slnsf 

A  .  33  and  I  are  given  in  Appendix  D. 
q  q  q 

An  important  cancellation  occurs  between  lim  J  and  the  first 


term  of  lim  J3 . 

€—>'0 

0 


€— »0 
o 


K2  .  No  singularity  is  involved  in  K2 .  Taking  the  limit 
value  before  integration  and  replacing  the  summation  with  in¬ 
tegration 


Ka  = 


00  c 

\  \ 

Jm=+l/2  J, 


a 

6 


dx 


5=0 


277 


a-iq0 


|dx2  ^X2+r2+p2-2rpcos0 


-iq0 


d0  dm 


X  =  m(  l/a)  9 


(38) 


0=0 


yx2+r2+p2-2rp 


•d© 


COS0 


X  =  2ieo 
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Using,  the  Lipschitz  integral6  and  Neumann's  addition  theorem7 


2tt 


G  =■ 


w 


-iq8 


d-0 


X2+ra+p2-2rpcos0 


(39) 


2rr 


0  =  0 


3-iq0 


,-|X|k 


k=0 


D  emJm(k:r)Jm(kp)cos  m0  dkdQ, 
m=0 


where 


e  = 
m 


1  for  m  =  0 
'2  for  m  7^  0 


Further, 


G  2^|  Jq(kr)  Jq(kp)e~  lXl  k  dk 
k=o 

=VF^Qq-l/2^Z^  ^ref- 


where 


=  r  +  p  + 

2rp 


( 41) 


The  function  Q  ,  /0  .  ,  .  -  „  ,  .  .  ,  ,  ,  , 

nq-1/2  is  a  spherical  function.  Numerical  tables 

of  Vl/2  necessary  for  the  following  treatment  are  given  in 


ref.  9. 


2  a  ,  ,  az 
3x  G  =  yF/if  az  ^q-1/2  z  ax 


2x  q-l/2 
(v^p)3  Z2-l 


z  Qq_l/2(z)  -  V3/2(z) 


(  >12) 
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The  behavior  of  ^y2  for  large  Z  is 


hi- 


n/o  =  1/2  X  f2(  n+s+1/2)  ^_z  ^n+s+1/2 

n-V^  s_0  i — : - -  '  7.4-n ' 


s  (n+s) 


Z+l' 


Hence  in  the  present  problem  with  q  >  0* 


ax 


G 


=  0 


u  j- 


(X=co 

Thus*  finally 

X  q-1/2 


K- 


(Vrp)3  Z^2-l 


Zo  %-l/2^Zo^ 


Qq-3/2^Zo^ 


where  „  =  o>‘ 

o 


2rp 


(43) 


(44) 


(45) 
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NUMERICAL  RESULTS  AND  DISCUSSION 

Calculations  have  been  performed  to  assess  the 
theoretical  approach  evolved  here  and  to  determine  the  three- 
dimensional  effects  in  the  stationary  and  nonstationary  flow 
cases  as  well  as  the  effects  of  such  important  parameters  as 
blade-area  ratio  and  pitch  angle »  The  calculations  have  been 
restricted  to  a  four-blade  propeller  with  sector  type  blade 
form  and  blade-area  ratios  A-g  =  0.22,  0.44  and  0.66  and  to  a 
constant  amplitude  gust  velocity  from  root  to  tip  expressed 
by  W  =  W  e_lq<^. 

The  results  are  exhibited  in  Figs.  4*  5  and  6  and 
summarized  in  Table  I  for  frequencies  q  =  0  to  4  and  for 
representative  values  of  inverse  pitch  angle,  a  ,  which 
bracket  the  range  of  pitch-diameter  ratio  of  practical  In¬ 
terest  j  i.e.j  a  =  2t r,  P/D  =  l/2  and  a  =  w,  P/D  =  1.  In  the 
same  figures  the  results  of  calculations  for  Weissinger's 
two-dimensional  model  by  a  stripwise  method  are  presented  for 
comparison. 

As  a  measure  of  the  three-dimensional  effect,  R3/2 
is  introduced  to  indicate  the  ratio  of  the  three-dimensional 
results  for  total  lift  to  the  corresponding  two-dimensional 
values 

R  _  load  by  three-dimensional  calculation 
3/2  load  by  stripwise  two-dimensional  approach 

Since  the  three-dimensional  and  two-dimensional  results  are 
obtained  through  the  same  Weissinger  approximation,  R3/2  can 
be  said  to  be  a  well-defined  measure  of  the  three-dimensional 
effect . 

It  is  seen  from  Figs.  4,  5  and  6  that  the  ratio 
R3/2  decreases  for  decreasing  order  of  harmonic  q.  This  in¬ 
dicates  that  the  over-all  correction  for  three-dimensional 
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effects  is  much  more  pronounced  in  the  stationary  case  than 
in  unsteady  flow,  a  fact  which  can  be  explained  at  least 
qualitatively,  in  the  following  manner.  The  propeller,  as 
any  other  lifting  surface  of  finite  aspect  ratio,  is  accom¬ 
panied  by  two  free  vortex,  systems.  As  a  result  of  the  span- 
wise  gradient  of  the  circulation,  a  free  vortex  system  is 
developed  along  the  helicoidal  surface.  On  the  other  hand, 
the  timewise  change  of  circulation  leads  to  a  radial  free 
vortex  system.  The  first  system  is  characteristic  only  of 
three-dimensional  flow  conditions  whereas  the  second  is  pre¬ 
sent  as  well  in  the  two-dimensional  model  and  to  the  same  ex¬ 
tent.  In  the  marine  propeller  case  due  to  high  load  and 
rather  low  pitch,  the  effect  of  the  first  vortex  system  is 
heavily  accumulated  in  stationary  flow  conditions,  thus  there 
is  a  large  discrepancy  between  two  and  three-dimensional  re¬ 
sults.  In  nonstationary  flow  conditions  the  helical  free 
vortex  changes  sign  q  times  per  pitch  and  hence  the  accumu¬ 
lated  effect  is  mitigated  considerably.  The  higher  the  fre¬ 
quency  the  smaller  the  effect  of  the  helical  free  vortices, 
therefore  the  closer  the  three-dimensional  results  to  the 
corresponding  two-dimensional  case. 

The  results  also  indicate  that  for  higher  blade- area 
ratio,  i.e.,  for  wider  blade  where  the  interaction  between 
blades  becomes  greater,  the  three-dimensional  effect  is  more 
pronounced.  In  fact.  Fig.  7  shows  that  for  larger  blade-area 
ratio,  coincidence  of  the  results  by  both  approaches,  three- 
dimensional  and  two-dimensional,  will  be  delayed  to  higher  re¬ 
duced  frequencies.  As  the  area  ratio  decreases,  coincidence 
will  occur  at  lower  reduced  frequencies;  in  other  words, 
three-dimensional  effects  become  less  important.  The  effect 
of  inverse  pitch  angle,  a  ,  is  not  as  important  (in  the  range 
■n  <  a  <  2tt)  as  is  the  blade-area  ratio.  However,  it  can  be 
said  that  the  smaller  a  the  closer  the  three-dimensional  load 
to  the  two-dimensional  value. 
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For  an  airscrew  propeller  with  small  blade-area 
ratio  ( small  solidity)  and  high  pitch  angle  ( small  a)  the 
three-dimensional  effects  should  be  less  important  than  for 
the  marine  propeller  with  larger  area  ratio  and  smaller  pitch 
angle.  This  conclusion  is  borne  out  also  by  the  comparison 
on  the  same  figures  of  the  phase  angles  between  downwash 
velocity  and  loading  obtained  by  both  approaches.  It  is  seen 
that  as  the  order  of  harmonics  increases,  the  three-dimension¬ 
al  value  for  the  phase  angle  tends  towards  the  corresponding 
value  for  the  two-dimensional  case.  For  smaller  blade- area 
ratio,  the  deviation  is  minimized  at  lower  reduced  frequencies. 

On  the  basis  of  the  present  investigation  revision 
can  now  be  made  of  the  method  used  in  ref.  3  to  obtain  vibra¬ 
tory  thrust  and  torque  and  the  large  discrepancy  shown 
there  between  theoretical  and  experimental  results  can  be  ex¬ 
plained.  The  method  of  ref.  3  was  based  on  the  work  of 
Ritger  and  Breslin2  who  used  two-dimensional  unsteady  airfoil 
theory  in  a  stripwise  fashion  in  conjunction  with  Burrill’s 
semi- empirical  correction  factors.11  Burrill’s  corrections, 
which  include  among  others  those  for  three-dimensional  effects, 
are  determined  for  stationary  flow  (  steady- state)  only.  Use 
of  these  corrections  in  refs.  2  and  3  for  both  steady  and  un¬ 
steady  flow  is  now  proved  invalid.  The  fact  that  the  calcu¬ 
lated  load  amplitude  in  unsteady  flow  was  found  to  be  l/3  of 
the  experimental  value  can  now  be  explained  by  the  results  of 
the  present  theory  which  show  that  the  correction  factor  for 
three-dimensional  effects  is  quite  different  for  stationary 
(q  =  0)  and  nonstationary  cases.  From  the  figures  and 
Table  I,  it  is  seen  that 

R3/2  for  q  =  0  _  l 

R3/2  for  q  =  4  ~  3 

which  indicates  that  the  values  of  the  unsteady  loads  ob¬ 
tained  by  the  method  of  ref.  3  were  overcorrected.  If  the 
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right  correction  factor  were  used  in  the  unsteady  two-dimen¬ 
sional  approach  the  theoretical  result  would  have  been  closer 
to  experiment. 


CONCLUSIONS 

The  unsteady  lifting  surface  integral  equation  is 
derived  for  the  marine  propeller  case  and  its  solution  ob¬ 
tained  by  using  the  Weissinger  approximation.  The  kernel 
function  is  expressed  in  closed  form  after  some  mathematical 
simplification.  The  applicability  of  the  Weissinger  method 
for  the  nonstationary  flow  case  is  studied  and  the  results 
indicate  the  validity  of  the  method  up  to  reduced  frequency 
1.3  which  corresponds  to  the  blade  frequency  of  a  four-bladed 
propeller  of  area  ratio  about  0.4  and  of  a  three-bladed  pro¬ 
peller  of  area  ratio  of  0.6. 

The  numerical  examples  restricted  to  a  four-blade 
propeller  with  sector  type  of  blade  form  of  area  ratio  0.22, 
0.44  and  0.66  and  to  a  constant  amplitude  gust  show  that  the 
difference  between  the  three-dimensional  approach  and  the 
two-dimensional  stripwise  approach  diminishes  as  the  order  of 
harmonic  increases,  as  the  area  ratio  of  the  blade  decreases 
and  with  increasing  pitch-diameter  ratio. 
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EXA  CT  INTEGRAL  PATH 
APPROXIMATED  INTEGRAL  PATH 


FOUR- BLADE  PROPELLER 
(N  =  4) 


FIGURE  2. 


APPROXIMATION  OF  INTEGRAL  PATH 


FIGURE  3.  DOMAINS  OF  INTEGRATION 
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FIGURE  4a.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  AND 
THREE-DIMENSIONAL  WEISSINGER  MODELS 


FIGURE  4b.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  AND 
THREE-DIMENSIONAL  WEISSINGER  MODELS 
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FIGURE  4c.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  AND 
THREE- DIMENSIONAL  WEISSINGER  MODELS 
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FIGURE  5a.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  AND 
THREE-DIMENSIONAL  WEISSINGER  MODELS 


FIGURE  5b.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO- DIMENSION AL  AND 
THREE-DIMENSIONAL  WEISSINGER  MODELS 


FIGURE  5c.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  AND 
THREE-DIMENSIONAL  WEISSINGER  MODELS 
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FIGURE  6a.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  AND 
THREE-DIMENSIONAL  WEISSINGER  MODELS 
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FIGURE  6b.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  AND 
THREE-DIMENSIONAL  WEISSINGER  MODELS 
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FIGURE  6c.  SPANWISE  LOAD  AND  PHASE  ANGLE  DISTRIBUTIONS  FOR  TWO-DIMENSIONAL  AND 
THREE-DIMENSIONAL  WEiSSINGER  MODELS 


FIGURE  7.  THREE-DIMENSIONAL  CORRECTION  FACTOR  R3/2  AND  PHASE-ANGLE  DIFFERENCE  BETWEEN 
THREE  AND  TWO-DIMENSIONAL  ANALOGS  AT  0.7  RADIUS  VERSUS  REDUCED  FREQUENCY  AT 
VARIOUS  BLADE  AREA  RATIOS  FOR  0  -  IT  AND  2.TT 
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APPENDIX  A. 

APPLICABILITY  OF  THE  WEISSINGER  METHOD 
TO  THE  NONSTATIONARY  PROBLEM 

In  order  to  clarify  the  applicability  of  the  Weiss- 
inger  method  to  nonstationary  problems.  Sears’  sinusoidal 
gust  problem  is  treated  by  this  method  in  this  section. 

Consider  the  system  shown  in  Fig.  A-l.  The  induced 
velocity  due  to  a  pulsating  dipole  located  at  (x',y’,z’)  is 
obtained  in  the  way  described  in  "Formulation  of  Problem, " 

( eqs .  1  to  8)  . 


W(x,y,z,t) 


elcot  a  a 

47T/iU  3z  dz  ’ 


x-x' 


i^{  T-x+x’) 


\Jtz+(  y~y ')  2+(  z-z  *) 


■  dr 


(A-l) 


It  must  be  noted  here  that  the  dipole  axis  is  taken  in  the 
negative  z'-direction  so  that  the  positive  direction  of  the 
lift  will  coincide  with  positive  z*.  If  the  dipoles  are  dis¬ 
tributed  uniformly  frpm  y’  =  -®  to  y'  =  on  the  1/4- chord 
line,  the  induced  velocity  upon  the  3/4-chord  line  is 


W(  t) 


lim 

ZO 


V 


loot 


,  +oo 


.  +1 


4TTjUU 


I  az2 

Y=-°° 


ik(  r-l) 


T  =  — oo 


-\Jt2+Y2+Z2 


dTdY  (  A- 2) 


where  Lq  is  the  lift  per  span  and  k  =  ^  is  the  reduced 
frequency.  The  integral  part  of  eq.  A- 2  is  written  as 


0 


lim 

Z->0 


-(-00 

-a2 

3Z2 


Y=-oo 


+1 


pikr _ 1 

\/t2+ Y2+Z2 


drdY 


lim 

Z-K) 


a 

3Z 


ikT_2Z_ 

tS+Z2 

—  03 


dr 


(A-3) 
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G  is  divided  into  three  parts :• 


G  =  2(G  +G  +G  ) 

12  3 


G  =  lim 
1  Z-»0 


G  =  lim 
2  Z-+0 


G  =  lim 
3  Z-»0 


a 

^Z 

r 

a 

3z 

f 

a 

oz 

r+1 

J  +£ 

ikr  Z  . 
e  -v — ~  dT 
r2+Z2 


elkT  -1-  dr 
t2+Z2 

ikr  Z 


r2+Z2 


dT 


(  A-4) 


where  it  is  assumed  that  |e|  «  1, 


Since  the  integration  range  of  G  includes  no 

/  a  \  l 

singularity  ,(  lim  )  can  be  taken  before  integration,  and  G 

\Z-»0  dZ/  : 

is  obtained  in  the  following  form: 

’~e  ikr  1 


G  = 
i 


_2 


dT 


1_ 

e 


cos  ke  ~kS  +  kSi(ke)  -i— sin  (ke)  +  lkCi( ke)  (  A- 5) 


where  Ci  and  Si  are  the  cosine  and  sine  integrals  which  are 
given  by 

cos  t 


Ci(x)  =  - 


x 


and  Si(  x)  = 


x 


sin  t 


dt 


dt  ,  respectively. 


G  = 

3 


In  a  similar  way. 


ikT 

dr 


R-9^0 
A- 2 


=  -  cos  k  +  cos  ke  -  kSi(k)  +  kSi(ke) 

-  i  sin  k  +  i  —  sin  ke  -  ikCi(ke)  +  ikCi(k)  . (A-6) 

In  the  Integral  region  of  G  ,  there  is  a  singular 

2 

point.  Therefore  the  limit  Z-^0  must  be  taken  last.  The 
exponential  is,  in  series  form, 

gikT  _  ^  +  _  i(kT)2  ....  (  A- 7) 

The  contribution  of  the  first  term  of  G  is 


G 

21 


lim 

2*0 


d 

BZ 


z 

t2+Z2 


dT 


(  A-8) 


That  of  the  second  term  (and,  in  general,  the  terms  of  odd 
powers  of  t)  is  zero,’  that  is. 


+e 

G  .  =  lim  1  ikr  — - —  dr  =  0 
22  Z-0  dZ  J  r2+Z2 

-e 


The  contribution  of  the  third  term  is 


G 

23 


lim 

Z-*0 


a 

Hz 


-  e 


Z 

t2+Z2 


dr 


(A- 9) 


(A- 10) 


Summing  up 
G  = 


=  -  k2e 

these  results  and  taking  the  limit  as  e--*0  leads  to 
lim  2  ( G  +G  +G  ) 

£-0  123 

2  [~eik  +  ik(  Ci( k)  +  iSi(k)  +  i  J)  ]  (A-ll) 
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which  upon  substitution  into  eq.  A- 2  leads  to: 


W(t)  =  eiuyk  [1  -  ike“lk  (Ci(k)  +  iSi(k)  +  1  J)  ]  (  A- 12) 

If  the  upwash  velocity  at  3/4~chord  point  is  given, 

the  complex  amplitude  of  the  lift  Lq  is  obtained  by  the  follow¬ 
ing  equation! 

L 

V  =__2_[1  -  ike“lk  (  Ci( k)  +  iSi(k)  +  i  J)  ]  (  A-13) 

3/4  27T/iU  * 

The  upwash  velocity  distribution  is  given  as  V(  x) e^0^  where 
V(x)  may  be  expanded  in  Birnbaum's  form 

00 

V(  x)  =  Vq  -  ^  ^  v^cos  na  (A-l4) 

n=l 

and  ^2/b  obtained  by  setting  x  =  1/2. 

However,  when  the  series  V( x)  is  composed  of  more  than 
two  terms,  the  straightforward  application  ^3/4  =  V(l/2)  leads 
to  very  poor  results  for  the  Weissinger  method.  This  diffi¬ 
culty  can  be  overcome  by  an  averaging  process!  is  taken 

as 


7T 

V3/4  =  M  V(x)  (1  "  cosa)  da  (A-15) 

^0  7T 

where  x  =  -  cos  a  and  the  integral  operator  ~  \  (l  -  cos  a) 
da  is  nothing  but  Glauert’s  lift  operator.  J  0 

If  eq.  A-15  is  used  in  the  stationary  state  case,  the 
Weissinger  approximation  gives  exactly  the  same  results  as 
thin  airfoil  theory.  On  expanding  V( x)  in  Birnbaum's  form 
and  applying  the  Glauert  integral  operator,  it  is  easily 
found  that  the  result  obtained  is  the  velocity  at  3/4-chord 
point  due  to  the  first  and  the  second  term  components  of  the 
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upwash  velocity  V(x)  in  the  Birnbaum  form. 

It  must  be  noted  once  more  that  the  application  of  the 
operator  gives  exact  results  only  in  the  s tationary  case. 
However*  in  the  following*  the  application  of  this  operator 
to  nonstationary  problems  is  tried.. 

In  Sears*  sinusoidal  gust  problem*  the  upwash  velocity 
V(  x*  t)  is  given  in  the  form  of 


V(  x*  t) 


and  with  x  =  -  cos  a 


( A-l 6) 


V(  x*  t)  =  V  eia)telk  cos  a 


(A- IT) 


where  k  =  ^  is  reduced  frequency.  If  the  integral  operator 
is  applied 


V3/4  =  Vo  * 


7T 

I* 

L 

7T 


(1  -  cos  a  )  elk  cos  a  da 


=  V 


7T 


l  ( 1  -  cos  a )  (  Jq(  k)  +  2  inJn(  k)  cos  n  ot )  d 

J  .  n=l 


a 


=  V  [J  (k)  -  iJ  (k)  ] 

0  0  1 

On  substituting  in  eq.  A-13 
Vo[Jo(k)  -  ±J1(k)]  -25PT 


(  A-l8) 


-  [1  -  ike"ik(Ci(k)  +  iSi(k)  +  iS)  ] 


(  A- 19) 

The  value  of  Lo/^7r/iUV'o  from  eq.  A-19  is  shown  in  Pig. 

A- 2.  In  the  figure*  the  exact  Sears’  value  is  also  plotted. 
The  agreement  between  them  is  fairly  satisfactory  in  the  re¬ 
duced  frequency  region  from  0  to  1.3.  The  value  of  reduced 
frequency  1.3  corresponds  to  the  blade  harmonic  ( q=N)  in  a 
propeller  of  area  ratio  about  0.4. 

R-940 
A- 5 


APPENDIX  B. 

GEOMETRICAL  MEANING  OF  THE 
SIMPLIFICATION  IN  THE  KERNEL 


The  simplified  expression  of  the  kernel,  eqs .  24, 
is  obtained  from  the  original  kernel  equation,  eq.  22,  by 
using  two  approximations,  one  in  the  integral  path  and  the 
other  in  the  direction  of  differentiation.  The  validity  of 
these  simplifications  becomes  clear  when  the  geometrical 
significance  of  eq.  22  is  compared  with  that  of  eqs.  24. 

Equation  22  indicates  the  following  manipulation: 
the  inverse  of  the  Descartes  distance. 


R[(  0,  r,  0)  ,(X,  p,  ©)]  \jxz+vz+pz-2vp  cos  0 


which  has  its  ends  on  a  group  of  helicoidal  surfaces 

(n  =  1,  2,  . ,N) ,  is  differentiated  at  both  ends  in  the 

directions  perpendicular  to  the  surfaces  and  then  integrated 
moving  the  end  (x,p,  0)  along  the  helical  trajectory  on  the 


helicoidal  surface  X  =  —  0 

a 


-  e  . 

a  n 


On  the  other  hand,  the  simplified  eq.  24  implies  a 
similar  procedure,  but,  this  time,  the  helicoidal  surfaces 
are  replaced  by  segments  of  planes.  Since  the  greatest  con¬ 
tribution  of  the  integral  in  eq.  22  and  eqs.  24  arises  in  the 
rather  narrow  region,  where  the  helicoidal  surface  in  eq.  22 
can  be  regarded  as  a 'plane,  eqs.  24  can  be  a.  very  good  ap¬ 
proximation  for  eq.  22,  because  in  this  region  the  relative 
direction  of  the  differentiation  with  respect  to  the  vector 
(0,r,0)  (X,p, 9)  is  the  same  in  both  expressions.  Numerical 
results  have  been  obtained  for  K2  when  r  =  .75  and  p  -  .85 

e0  30o 

and  X  is  taken  from  ^  to  -g-g-,  derived  from  both  expressions 
and  from  eq.  38  where  the  summation  is  replaced  by  integra¬ 
tion.  They  are  compared  in  the  following: 
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K2  for  q 

=  0,  r  = 

.75,  p 

=  .85 

a 

=  2n 

a 

=  7T 

Eq. 

22  from  X  = 

eo 

2a 

to 

30o 

2a 

- 

25.1 

- 

9.6 

Eqs 

.  24,  m  =  1 

- 

ro 

- 

9-0 

Eq. 

38  from  X  = 

0o 

2a 

to 

360 

2a 

- 

23.8 

- 

10.8 

The  agreement  among  them  is  satisfactory  for  practical  pur¬ 
poses  . 

The  values  of  from  eq.  22  and  from  eq.  25  de¬ 
rived  from  24  are  compared  below  for  the  same  values  of 
r  =  .75  and  p  =  .85.  (The  half-angular  chord  -  ®0=  10° 

for  the  calculations.) 


q 

Kjeq.  25) 

a  = 

Kx( eq.  22) 

2ir  a  =  3tt 

2 

a  =  tt 

0 

228.0 

218 

23.4 

205 

1 

224.3  -  5.7i 

216 

-  6i 

4 

200.8  -  19. Oi 

194 

-  l8i 
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APPENDIX  C. 

EVALUATION  OP  INTEGRAL 


ih 

i*  >  r 


U 


>±2qf 


*1 


(\}l-k2sin2^/)3 


df 


Watkins10  derived  the  recurrence  formula  for  the 
special  integral  limits  f  =  S  and  f.  -  0.  A  somewhat 

U.  £—  <V 

similar  method  gives  the  formula  for  general  integral  limits, 


In  the  first  place*  the  following  notations  are  in¬ 


troduced: 


v,n 


& 


V, 


ei2n t/j 


( l-k2sin2^) 


W  a* 


'  K,n  +  1  13 


u 


vs  n 
cos2n^ 


■t 


% 

The  expression 


( l-k2sin2^) 


d^  +  i 


(0-1) 


sin2'n^ 

— :"Ivv/2  a* 


t 


( l-^sln2^) 


l 


cos  2(n+l)^/+cos  2(n-l)^=cos  2n^(2-^-)+-^  cos  2n^(  l-^sin2^) , 

(C-2) 

(C-3) 


k£ 


gives  the  relation 


-c  +  Ic 
‘  a  ,  n+1  3,  n-1 


ir  +  r.r  _  ,=  2(1- A)  i°  ^  +  jl  i 

k2 


4  Tc 


n  k2  i?n 


Also,  if  the  relation 


sinicos^ 


.  d^  =  — 


e 

J  ( l-k2sin2^)  k2  ( l-k2sin2^) 

is  used*  the  following  expression  is  obtained: 


(c-4) 


R-9^0 

C-l 


-It  ,  =  -4 


3,n+l  s,n-l 


sin2n^ 


_k2  (l-^sin2^)1/2 


ty. 


u 


2n  Tc 

?  1 


t 


(0-5) 


l 


Eliminating  I,  from  eqs .  C-3  and  C-5, 

l9n 

( 2n-l)  1°  ^  +  (2n+l)  Ic 


a,n+l  1  ^3,n-l  v2'  “3,n 


=  4n(l  --L)  1° 
k2 


sin2n ip 


^  ( l-k2sin2?,i/)  4/2 


t 


u 


(C-6) 


l 


In  a  similar  way., 

r  s 


<2n-D  T-l.n+l  +  <2n+l)  Is,n-1  “ted-  p  li>n  (0-7) 


J±_ 

k2 


cos2n^ 


(.  l-k2sin2^) 


l/2 


Thus,  considering 


I  =  I 

3iQ  q 


( C-8) 


and  using  the  relation  of  eq.  C-l,  the  recurrence  formula  for 
as 

te_  ( 1  -  _L)  1  +  Sail 


I  is  written  as 

q 


q+1 


2qTl 


k" 


q  2q-l  q-1 


-  i 


,12q^ 


2q_1  k2  (i-k2sin2^)l/2 


t 


(C-9) 
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In  order  to  use  this  formula,  the  evaluations  of  I„  and  I 

J  0 

are  necessary.  Some  elementary  calculations  give 


1 

1-k2 


E(k^) 


k2  U/2)sin2^ 

( l-k2sin2^) 


*1 


P(k,^) 


t 


and 


_2_  1 _ 

k2  ( l-k2sin2^) 1//2 


t 


u 


i 1 


t 


(C-10) 


where  F(k,^)  and  E(k,^/)  are  the  first  and  the  second  kind  of 
incomplete  elliptic  integral.  Combining  the  cosine  and  the 
sine  terms. 


I  =  I 
■L0  x3,0 


1-k2 


e(  k,^)  -  ks 


(l/2)  sin2 f 


4l-l 


k2sin2^J 


Vi 


(C-ll) 


I  =  I 

II  13,1 


=  (1  -  JL>  ln  +_  P(k ,ip) 


kc 


f. 


u 


+  i  — 


^  ( l-k2sin2^) 


(C-12) 
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APPENDIX  D. 
INTEGRATION 


The  contribution  of  the  small  area  including  the 
singular  point  0=0,  p  =  r  is 

J u 


=  lim 
x-*0 


V€o 


P~rp~eO J  0=  -7 


giqe 


clx2  y^x2+r2+p2-2rp  cos© 


d0  dp 


(D-l) 


As  is  stated  in  "Evaluation  of  The  Kernel  (K  )"  on  p.  10, 

the  case  of  -y|«  1  ,  e  «  r  is  considered.  Therefore 
•  1  o  p 


=  lim 
x-»0 


'0  1+7 


dx2  Vx2+e2+r>  2  ©2 
e0  J-y  P 


de  d0 


(D-2) 


and  taking  the  lim  finally 
x->0 


J,  =  -  4  — 
1  r 


1 

\<v 


PHliyy)2  ( eQ) 2 


(D-3) 


The ref ore. 


lim  Ji  =  "  4  FT 


£T>0 


P  0 


(D-4) 
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APPENDIX  E.. 
INTEGRATION  J . 


The  integration  domain  does  not  include  the  singular 
point,  therefore  can  be  evaluated  by  the  recurrence  re¬ 
lation  obtained  in  Appendix  C.  However,  the  limiting  process 
eQ— >-0  precludes  using  a  numerical  method  for  the  p  integra¬ 
tion.  In  this  section  this  difficulty  is  overcome  by  evalua¬ 
ting  the  integral  in  an  approximate  but  analytical  form  made 
possible  because  /3X,  j8a  «  r  ■ 


The  limiting  process  is  taken  first,  then  In  ellip¬ 
tic  integral  form,  _ 


Ja=( -l)q+1  2 


Ve0rrp+P2 


+ 


P+£0J 


1  -  0n 

4(7T+en+-^) 


1  N 

(  rp+P) 3  £ 


eig(^  dfdp 


n=l|  \j  l-k2sin2^ 


4r  p 

where,  k2  =  y—J— 

(rp+p)2  • 

Further,  J  is  divided  into  two  parts. 


(E-l) 


J31=(-l)q+1  2 


J32=(-l)q+1  2 


>rp  E0 

r  +/3~  -j 

p 

+ 

L-V't  ' 

r  +eA_ 
P  o 

rfve°  , 

.rp+fe- 

+ 

-8, 

P 

rp+€0 

—7 r 


si2q  f 


(r  +p)3j  J(l-k2sin2^): 

y  0 


dfdp 


( E-2) 


(  ^p+p) 3 


-  ?(7r-0b)  -7T 


“  2"(  TT-t-g") 


(eq.  E-2  con’t.  on  next  page) 
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E-l 


e 


N 

+  £ 

n~2 


0\n 


i(7r+Veb)J 


^12qf 


^l-k2sin2^/ 


•  d^dp 


(E-3) 


with  k2  =  and 

( ^p+p) 2 

ks  =-yJl  -  k2 
J31  can  be  written  as 

Jai=  Jai  +  Jat 


fix 


J3l  =  (-1) 


q+1 


2  r2v^ 


(2rp)« 


( 1+k ' )  I»  dk' 


2Veo 


Jat  *  (-l)q+1 


P 


2rp+eo 


(E-4) 


-7T 


,i2q^ 


where  I*  =  I  , - c 

^  J  q  l-k2sin2^ 


-df 


(E-5) 


The  recurrence  relations  obtained  in  Appendix  C  can  be  adapted 
for  the  elliptic  integral  I*  .  Then 
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E-2 


q+1 


4 


q  i  1+k'2 


2q-l  4l-k' 


■)  + 


2q+l  p 
2q-l  x  q-1 


I'"  -  -[■ 


0 


k 


fa  E 


(vi-k’2)] 


(e-6) 

(e-7) 

(E-8) 


where  K(  k)  and  E(  k)  are  complete  elliptic  integrals  of  the 
first  and  second  kind.  These  formulae  give 


1 1  =  - 
x  0 

;[^E] 

z\  - 

« k%  * 

- ^75-  K 

1-k'2 

] 

z,2  =  - 

■  [(^- 

3)  i^E- 

49  r 

4 

-k'2 

I?3  " 

'<?  «•  ■ 

-f «)  -4 

_  E  - 

5 

(32 

If4=  - 

"(i§§Q‘- 

i44^) 

22e 

k'2 

-  0 

(E-9) 


-  4)  — ^  K 

5  1-k’2 


1+k'2 
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In  the  following  treatment,  since  j31,  <  <  r^  and 

hence  k*  <  <  1,  the  terms  smaller  than  the  zero  power  of  k* 
in  the  integrand  of  eq.  E-4  are  neglected.  The  elliptic  in¬ 
tegral  E  and  K  can  be  expanded  in  powers  of  k'2  as 

E  -  1  +  |  (to  h  -  |)  k's  +  it  -  +  .  .  . 

K  *  W  +  ¥(to  F  -  1}  k'‘  +  Uln  £  '  6)k'4  +  '  •  ’ 

(e-io) 

Hence,  from  eqs.  E-9?  a  sufficient  form  of  1'^  for  the  above- 
mentioned  approximation  is 
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Eq.  E-4  becomes 
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If  /3  and  /3a  are  chosen  to  satisfy  the  relation 
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lim  J  is  given  in  the  following  form: 
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As  is  seen  in  eq.  E-3,  even  the  point  in  the  inte¬ 
gration  domain  of  J32  which  is  closest  to  the  singular  point 
is  at  a  distance  from  it  equal  to  r  0^,  where  0^  is  the 
angular  semi-chord. 

If  /3X  and  02  are  very  small  compared  to  r  0^,  k 
can  be  assumed  unity  throughout  the  whole  integration  domain 
of  J0.  Thus, 

32  ’ 
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The  recurrence  relation  obtained  in  Appendix  C  holds  also  for 

4  e12^ 


I'1 'q,  but  in  this  case  it  becomes 
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where  the  double  sign  should  be  taken  in  the  following  senses 
+  for  the  region  cos  f > 0 
-  for  the  region  cos  f<  0 
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FIGURE  A-2.  APPLICATION  OF  WEISSINGER  METHOD  TO  SEAR'S  GUST  PROBLEM 
FOR  NON -STATIONARY  CASE 
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